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We examine the relevance of magneto-elastic coupling to describe the complex magnetic and 
structural behaviour of the different classes of the iron superconductors. We model the system as 
a two-dimensional metal whose magnetic excitations interact with the distortions of the underlying 
square lattice. Going beyond mean field we find that quantum fiuctuation effects can explain two 
unusual features of these materials that have attracted considerable attention. First, why iron 
telluride orders magnetically at a non-nesting wave- vector (7r/2,7r/2) and not at the nesting wave- 
vector (tt, 0) as in the iron arsenides, even though the nominal band structures of both these systems 
are similar. And second, why the (7r,0) magnetic transition in the iron arsenides is often preceded 
by an orthorhombic structural transition. These are robust properties of the model, independent of 
microscopic details, and they emphasize the importance of the magneto-elastic interaction. 

PACS numbers: 74.70.Xa, 74.90.+n, 75.80.+q 
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Introduction. — The recently discovered iron supercon- 
ductors with unusually high transition temperatures ex- 
hibit a rich phase diagram that includes structural, mag- 
netic and superconducting transitions [HQ- As such 
these materials are the latest playgrounds to study how 
in complex materials different phases compete, and how 
this unconventional setting eventually gives rise to super- 
conductivity. Theoretically, one of the goals at present 
is to identify the microscopic interactions that give rise 
to the rich phase diagram. This motivates us to present 
a microscopic study of the simplest symmetry-allowed 
model Hamiltonian describing magneto-elastic interac- 
tion. The results allow us to argue that this coupling 
contains physics relevant for the iron superconductors, 
and therefore it is an important microscopic ingredient. 

Crystallographically these materials have a layered 
structure, which is reflected in their energy bands with 
weak dispersion along the c-axis compared to that along 
the a6-plane [1, 01 ■ Consequently, it is often simpler to 
consider them as two-dimensional systems weakly cou- 
pled along the c-axis. The undoped and the lightly doped 
compounds usually undergo magneto-structural transi- 
tions from paramagnetic metals with tetragonal crys- 
talline symmetry to low temperature antiferromagnetic 
(AF) metals with either orthorhombic (in case of the 
FeAs systems) or monoclinic (in case of Fci+^Te) struc- 
tures. These transitions are suppressed in favour of su- 
perconductivity when they are either doped or put under 
external pressure. The AF order of the FeAs systems is 
at the wave- vector (tt, 0) in the Brillouin zone defined by 
the plane of the Fe atoms with IFe/cell, and it is often 
preceded in temperature by the structural transition. On 
the other hand, the AF order of Fei+yTe is at (7r/2, 7r/2), 
and the lattice distorts simultaneously. 

While it is likely that the (tt, 0) order of FeAs is a 
consequence of the Fermi surface nesting in these multi- 
band systems 0] , from the perspective of a band picture 



there remains at least two important puzzles concerning 
the magneto-structural properties. First, why Fei+yTe, 
whose nominal band structure is similar to that of the 
FeAs systems orders at (7r/2, 7r/2) and not at the nest- 
ing wave- vector (tt, 0). And second, why the AF transi- 
tion of the FeAs systems is often preceded in temperature 
by a tetragonal-orthorhombic structural transition. The 
main results of this study are to show that both these fea- 
tures are natural consequences of quantum fluctuations 
induced by the magneto-elastic interaction. We find that 
the former property is due to spin fiuctuations scattering 
with short wavelength phonons, and the latter is driven 
by critical spin fiuctuations near a (tt, 0) AF transition. 

The above questions, as well as the general magneto- 
structural properties of these materials, have been ad- 
dressed earlier from various points of view. Some of these 
are based on itinerant electron models which emphasize 
the physics of nested Fermi surfaces 0] . There are also 
studies that suggest that the electron-electron interac- 
tion is strong [7|, which justifies describing the magnetic 
properties by Heisenberg spin models with appropriate 
couplings 0. The structural transition has been viewed 
as a consequence of various kinds of electron order such 
as orbital ordering [9^ , spin nematic ordering fio t at tem- 
peratures above the magnetic transition ll|, Il2i | , as well 
as the ordering of orbital currents [isj . 

In the past there has been few studies of the magneto- 
elastic properties of the FeAs materials which concen- 
trated on the c-axis motion of the As atoms, and its 
strong influence on the Fe-As bond and eventually on the 
magnetism 1J| . These are motivated by the observation 
of a pressure-driven first order volume collapse transition 
in CaFe2As2, which is concomitant with the loss of mag- 
netism [l^. On the other hand, a microscopic study of 
the influence of the a6-plane distortions of the lattice on 
the magnetic sector using an effective model is currently 
lacking. This is despite the fact that the in-plane distor- 
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FIG. 1; Magneto-elastic coupling [Eq. ([T])] describing scat- 
tering of paramagnons (straight line) (a) with phonons (wavy 
line) at finite k, and (b) with orthorhombic distortion (dash 
line) at k ^ limit. The associated matrix elements gc,(k,q) 
and /i(q) are defined in Eqs. @ and ([HJ respectively. 



tions have an important symmetry-allowed coupling with 
the order parameters of the magnetic transition, which 
can make the transition weakly first order [l^ [I^ . 

Here we perform such a microscopic study using a 
model Hamiltonian Hme introduced phenomenologically 
in Ref. Il8l . and where it was argued that a mean field 
treatment of T-Lme describes several salient features of 
the magneto-structural transitions in the FeAs [l^ and 
Fei+j,Te 18| systems. In the case of FeAs it explains 
(i) why the magnetic and structural transitions are sep- 
arate in some materials while they are concomitant in 
others, (ii) why the transitions appear to be first order 
in some cases even though they are allowed to be sec- 
ond order from a symmetry point of view, and (iii) why 
the systems prefer a coUinear magnetic state instead of 
a non-collinear order [l^. In the case of Fei+yTe it ex- 
plains why in the magnetic phase the system undergoes a 
uniform monoclinic distortion, and a modulated one with 
wave-vector (7r,7r) such that the Fe-Fe bonds are alter- 
nately elongated and shortened (3). Our aim here is to 
go beyond mean field, and to study the quantum fluctu- 
ations of the elastic and the magnetic degrees of freedom 
and thereby examine further the relevance of Hm e ■ 

Model. — We consider a two-dimensional metal on a 
square lattice with magnetic and elastic degrees of free- 
dom that are coupled by the Hamiltonian 
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Here Aq is the coupling constant with dimension of energy 
per length, = u{ri) is the displacement of the Fe atom 
from its equilibrium position at r^, (ij) implies nearest 
neighbour sites, hij is the unit vector along the j— j bond, 
and Mi is the local magnetization at r^. We neglect spin- 
orbit coupling, and magneto-elastic terms of 0{u^) and 
higher. In the context of insulating magnets Hme de- 
scribes the variation dJ/dr of the Heisenberg exchange 
J with bond length, and its effects are well-studied [20j. 
In the case of metals, where such couplings are far less 
studied, Hme describes the bond- length dependence of 
the parameter Vij that characterize the nearest neigh- 



bour interaction VijPi-i^pji, with pia- being the density of 
electrons with spin a at site i. 

We study the system from the paramagnetic side, 
and we approximate to be 0(3) variables describ- 
ing spin fluctuations (paramagnons). Conceptually, they 
can be introduced as Hubbard-Stratonovich fields to de- 
couple the appropriate electron-electron interaction, af- 
ter which the electrons can be integrated out 2l|. The 
resulting action for the magnetic sector can be written 
as Sm = Eq,i-„,aX(7^(q>i'n)M*(q,z/„)Ma(q,:^„), where 
a = {x, y, z) in 0(3) space, M(q) is the Fourier trans- 
form of M(ri), and i/„ is a bosonic Matsubara frequency. 
We assume that the paramagnon propagator describes 
damped dynamics with Xo^i^^'^n) — no(q) -I- \vn \ , and 
having the bare dispersion Ho(q) — flo+^i [cos qx cos Qy — 
{cos2qx + cos2qy) /8] with f^o > i^i > 0. With this choice 
Ho(q) has a global minima at q = (tt, 0) (and at the sym- 
metry related points) which models the nesting property 
of the underlying Fermi surface. The approach towards 
the magnetic instability can be described by lowering fio- 
As our results are independent of the microscopic de- 
tails of the elastic sector, it is sufficient to describe it us- 
ing the simplest model compatible with the square sym- 
metry. The lattice variables are defined by the strain ten- 
sor eafjivi) = Uap + i/2 J2k=io[kaUf3{'k) -t- a o /3] exp(ik ■ 
Ti), where a,f3 — {x,y) in two dimensions. This includes 
the uniform strains Mq^, and the phonons described by 
u(k) which is the Fourier transform of u(ri) [13] • The en- 
ergy per unit area due to is given by Es = Cii{uxx + 
UyyY /2 -\- Ci2UxxUyy + Cgeu^j^/S, whcrc the Cs are the 
elastic constants in Voigt notation. The order param- 
eter for the tetragonal-orthorhombic transition is Uo = 
Uxx — Uyy, and the associated elastic constant is Co = 
(C\\ — C\2) 12. Next, the phonons are described by the ac- 
tion Sp = X;k,i.„,^t^/7Hk, i^n)w^(k, z^n)w^(k, i^„), where 
p is the polarization index. We assume the propagator to 
have the standard form L)~^(k, Vn) 



with L0\^^ the phonon dispersion. We obtain it from a har- 
monic theory in which the eigenvalues of the dynamical 
matrix Nap{\i) give mw^^, m being the atomic mass of 
Fe. We take Nxx = 4Cii sin^(fc:r/2) + Cee sin^(fcy/2). 



Nxx{x o y), and 



N„ 



(4Ci 



Cee) sin(fca;/2) sin(fcj^/2). Thus, the elastic sector is en- 
tirely characterized by the Cs and m. 

Before performing the calculations it is instructive to 
re- write the interaction in Fourier space as 

nME=2i\o 3a(k,q)u„(k)AC(k + q)Afa(q), 

a,a,k,q 

where ga (k, q) are matrix elements with 

gaiKi) ^ sm{ka + Qa) - siniqa)- (2) 

They play an important role due to the property that 
(/ct as fc 0, which implies that the coupling to 
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FIG. 2: Magneto-elastic coupling induced quantum fluctua- 
tions producing (a) static paramagnon self-energy 112 (q) due 
to scattering with phonons, and (b) correction to the or- 
thorhombic elastic constant Co due to paramagnons. 



the phonons vanishes in the limit of uniform displace- 
ments. In this limit the paramagnons couple to the clas- 
sical variables Uq/j, and in particular the coupling to Uq 
is associated with the matrix element 




(0,0) (71,0) (it.Tt) (5,5) (0,0) (0,0) (7t,0) (7C,IC) (5,5) (0,0) 



/i(q) = cosqx — cosqy. (3) 

The coupling of the paramagnons to the phonons (k ^ 0) 
and to Uo (k — > 0) along with their respective matrix 
elements are shown graphically in Fig. [T] 

As our model is two dimensional, strictly speaking it 
cannot be used to study the finite temperature magnetic 
transition due to the Mermin- Wagner theorem. In the 
following our strategy is to perform calculations at zero 
temperature (T) where the theorem is inapplicable, and 
to infer finite-T consequences using adiabaticity argu- 
ment. In general, at finite-T the magnitude of the results 
below are enhanced due to the thermal fiuctuations. 

Results. — (i) First we examine the static paramagnon 
self-energy obtained due to virtual scattering with the 
phonons. This is given by (see Fig. [2^) 

n2(q) = -4A^T (2™c.kM)"'5a(k,q)5/3(k,q) 

X eM(k)e(^(k)i:';,(k, ij^„)xo(k + q, ii^n), 

where £^,0^) are the polarization vectors. We perform 
the k-sum numerically, taking (Jloj^^i) = (10,2) meV 
(the precise value of these parameters do not affect the 
result qualitatively). We also use (Cn, C12, Cee) — 
(0.484,0.161,0.141) meV/A^ f^Qj^ elastic constant 
measurements on BaFe2As2 j22l | suitably normalized for 
two dimensions. In Fig.[3K we plot 112 (q) along high sym- 
metry directions of the Brillouin zone for various values 
of the coupling Aq. Note that the magnitude of 112 (q) 
increases with Aq as expected, but more importantly, 
it is peaked at {tt/2,tt/2) for all Aq (below we explain 
why). The consequence of this property is demonstrated 
in Fig. ^jp where we show the corresponding renormal- 
ized paramagnon dispersions (second-from-top to bottom 
curves) given by X~^(q) = no(q) ~ 112 (q), as well as the 
bare dispersion no(q) for comparison (topmost curve). 
Since n2(q) is always peaked at (7r/2,7r/2), for Aq > A^ 
(e.g., the bottommost curve in Fig.|3)3) the global minima 
of X~^(q) changes from the nesting driven wave- vector 



FIG. 3: (colour online), (a) Paramagnon self-energy 112 (q) 
along the high symmetry directions of the Brillouin zone for 
increasing values of magneto-elastic coupling Aq = I,-- - ,5 
(bottom to top) in unit of 10 meV/A. It is always peaked 
at (7r/2,7r/2). (b) The corresponding renormalized param- 
agnon dispersions X~^(ci) = no(q) ~ n2(q) (second-from- 
top to bottom curves respectively), and the bare dispersion 
Ilo(q) (topmost curve, shifted up by 0.2 meV for clarity). At 
large enough Aq (bottom-most) the global minima changes to 
(7r/2,7r/2) (dash arrow) from the nesting driven (vr, 0) wave- 
vector (solid arrow). 



(tt, 0) to the magneto-elastic coupling driven wave- vector 
(7r/2,7r/2). This implies that in this model the mag- 
netic instability is at (7r/2,7r/2) for sufficiently large Aq. 
Thus, the (7r/2, 7r/2) magnetic order observed in Fei+^Te 
can be explained as a signature of strong magneto-elastic 
coupling, while in the FeAs systems this effect is presum- 
ably weaker and the nesting driven (tt, 0) instability is 
preferred. Indeed, several authors have argued that the 
electron-electron interaction effects, which favours the 
magneto-elastic coupling, is stronger in Fei+j,Te than in 
the FeAs systems . 

The above result can be understood simply from the 
following argument. Since ga vanishes for /c — >■ 0, the k- 
sum above is dominated by large k (short wavelength 
phonons), and typically Wk/i ^ Eq the phonon band- 
width, which is the largest energy scale in the model. 
Thus, TJ2^^D^(k,iVn)xoO^ + (l,ivn) ~ 1/i^o, and its 
q-dependence is negligible. Therefore, the dominant q- 
dependence of 112 (q) is from the matrix elements gadp- 
Then a simple power counting argument, and the fact 
that in a stable lattice the diagonal components of NajS 
are larger than the off-diagonal ones, establishes that 

n2(q) = A + B (sin^ + sin^ Qy) — C cos cos Qy, (4) 

with A, B > C > 0. Note that this result is independent 
of the details of the paramagnon and the phonon dis- 
persions, and its main ingredient is the fact that Hme 
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involves the nearest neighbour sites. Including couplings 
with longer range will not change the result provided they 
are smaller in magnitude compared to Aq. This is plau- 
sible, since the electron-electron interaction is expected 
to reduce with distance. However, it needs to be checked 
using first principles calculation. Note also that the value 
of Aq depends intricately on the microscopic details. In 
particular, it depends on the relative strength of the nest- 
ing driven minimum at (7r,0), and on the magnitude of 
the longer range magneto-elastic couplings. 

(ii) Next we compute the correction to the orthorhom- 
bic elastic constant Co due to quantum fluctuations of 
the paramagnons which is given by (see Fig. [2}d) 

5Co = -i\lTY,h\ci)xl{ci,iVr,). 

For two dimensional spin fluctuations in the vicinity of a 
(tt, 0) magnetic instability (which is tuned by J7ci) we get 

<5Co = -4A^/(7r2l]i)log(l/(5), (5) 

where 5 — (ilo — ilciO/f^i measures the closeness to the 
magnetic instability at fio = i^cr- This implies that, 
before the system becomes magnetic (J7o > ilcr), the 
renormalized elastic constant (Co -I- 5Co) vanishes and 
the system encounters an orthorhombic instability. This 
softening of the lattice with SCq oc log(l/xm), where 
Xm is the dimensionless static magnetic susceptibility at 
(tt, 0), is in principle verifiable from neutron scattering 
data and from elastic constant measurements. A similar 
softening of Cq due to nematic fluctuations is discussed 



in Ref. |12|. Using adiabaticity argument it is possible 



to extend the phase boundaries to finite temperatures in 
the T — plane, and the resulting phase diagram shows 
that the magnetic transition is preceded in temperature 
by the orthorhombic instability. However, close enough 
to the finite-T magnetic transition the critical theory is 
non-Gaussian, and the quantitative aspects of the current 
treatment becomes invalid. 

In practice, even at T = and in three dimensions the 
log divergence is cutoff by the paramagnon dispersion 
along the c-axis, and the precise fate of the structural 
instability depends on microscopic details. Nevertheless, 
for quasi-two dimensional paramagnons one should ex- 
pect considerable orthorhombic softening, and therefore 
this mechanism explains why often (but not always) the 
(tt, 0) magnetic transition in the FeAs systems is accom- 
panied by an orthorhombic transition. In passing we note 
that a similar monoclinic softening is expected near a 
(7r/2,7r/2) magnetic transition (and can be relevant for 
Fei+yTe), but to capture this physics one needs to gener- 
alize I-Lme and include next-nearest neighbour couplings. 

Conclusion. — We studied the effects of quantum fluc- 
tuations induced by magneto-elastic coupling [Eq. ([T])] in 
a two-dimensional metal on a square lattice. The cou- 
pling describes the simplest symmetry-allowed scatter- 
ing between the paramagnons and the distortions of the 



lattice. At a qualitative level the model explains (i) why 
Fei+j,Te orders magnetically at (7r/2, 7r/2) and not at the 
nesting wave- vector (tt, 0), and (ii) why the FeAs systems 
often undergo an orthorhombic transition in the vicin- 
ity of the (tt, 0) magnetic transition. The former result 
is due to paramagnons scattering with short wavelength 
phonons, and the latter is driven by the critical spin fluc- 
tuations. We hope these results will stimulate further 
studies of the coupling using a more realistic model and 
in conjunction with first principles calculations. Such 
couplings can give rise to qualitatively new physics, and 
they can be relevant for correlated metals in general. 

The author is very thankful to E. Boulat, A. Cano, 
H. Capellmann, E. Kats, I. Vekhter, and T. Ziman for 
insightful discussions. 
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